Anomalous fluctuations of two-dimensional Bose-Einstein condensates 
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We investigate the particle-number fluctuations due to the collective excitations created in a two- 
dimensional (2D) and a quasi-2D Bose-Einstein condensates (BECs) at low temperature. We find 
^ that the fluctuations display an anomalous behavior, i. e. for the 2D BEC they are proportional to 

N 2 , where N is the total number of particles. For the quasi-2D BEC, the particle-number fluctuations 
are proportional not only to N 2 but also to the square root of the trapping frequency in the strongly- 
<• confined direction. 
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The remarkable experimental realization of Bose-Einstein condensates (BECs) has stimulated intensive theoretical 
3 i and experimental studies on weakly interacting Bose gases [1] . Much attention has been paid to the investigations of 
the thermodynamic properties especially the particle number fluctuations of three-dimensional (3D) interacting Bose 
gases [2-7]. Recently, quasi-lD and quasi-2D BECs have been realized [8,9], which with no doubt provides many 
new opportunities to explore the fascinating quantum statistical property of macroscopic quantum systems in low 
dimensions. For a 2D Bose gas, the particle number fluctuations due to the exchange of atoms between condensate 
O ] and thermal atoms have been addressed recently in [10,11]. 

It is well known that for the temperature far below the critical temperature, the collective excitations play a 
\ dominant role for the fluctuations of BEC. The fluctuations of condensate due to collective excitations in 3D Bose 
■ gases have been discussed in [2,5,12]. The scaling behavior of a 2D interacting condensate confined in a box was 
investigated in [5]. In the present paper, we shall discuss the fluctuations originated from collective excitations for 
2D and quasi-2D interacting condensates confined in a magnetic trap. 

For the Bose-condensed gas confined in a magnetic trap, the total number of particles N in the system is con- 
served and hence a canonical (or microcanonical) ensemble should be used. Within the canonical ensemble a fairly 
fSj [ general method has been developed recently [12,13] for studying the thermodynamic properties of the interacting 
Bose-condensed gases based on the calculation of probability distribution functions. The purpose of this work is to 
investigate the particle number fluctuations of a 2D BEC due to the collective excitations within a canonical ensemble. 

For this aim we first extend our previous theory developed in [12,13] on particle number fluctuations in BECs by 
including an effect of the quantum depletion. For a Bose gas confined in a trap, based on Bogoliubov theory [14] 
one can obtain that, at low temperature, the total number of particles out of the condensate due to the collective 
£^ ' excitations is given by 
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where (Nt) and (No) are respectively the number of excited atoms and atoms in the condensate. f n i — 
1/ (exp (/3e n i) — 1) is the average number of collective excitation nl present in the system at thermal equilibrium 
(f3 = l/(fcsT), fcs is Boltzmann constant and T is temperature) , while e n i = fuWnl is the energy of the collective 
mode characterized by the quantum numbers n and I. In the above equation, the term J v^dV represents the effect 
of the quantum depletion which does not vanish even at T = 0. The quantities u n i and v n i are determined by the 
following coupled equations: 
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where V ex t (r) is the trapping potential confining the Bose gas, n and g are respectively the chemical potential of the 
system and inter-atomic interaction constant, n (r) and no (r) are the density distribution of the Bose gas and the 
condensate, respectively. From Eq. (1) we see that (N n i) can be taken as the average number of the excited atoms 
related to the collective mode nl, while (N^i) — fnl is the average number of the collective excitation nl. From the 
form of Eq. (1), we see that the ratio between (N n i) and (N%) is given by 

_ iM - (f<i dV + f V ll dV )fnl+fv 2 nl dV 

Within the canonical ensemble, the partition function of the system with N atoms takes the form 



Z[N] exp 



(5) 



In the above equation, we have omitted the interaction between collective excitations. The prime in the summation 
represents the condition that the total number of atoms in the system should be conserved within the canonical 
ensemble, i. e., X)ni#o N n l — Snj^o VniN^i = N — No. One should note that in Eq. (5), is the occupation number 
of the collective excitation nl and should be integer number. For the convenience of calculations, by separating out 
the ground-state nl = from the state nl ^ 0, we have 

N 

Z \- N \= Yl i ex P l-PNoeo] Z (N T )} , (6) 
N =a 

where Zq (Nt) stands for the partition function of a fictitious system comprising Nt — N — N excited atoms which 
takes the form: 

Z (N T )= Yl exp[-t3N*e nl \. (7) 

The free energy of the fictitious system is A (Nt) = — fc^Tln Z (Nt). 

By using the developed saddle-point method proposed in [12,13], we introduce a generating function Go (T, z) which 
is given by 

oo 

G Q (T,z)= ]T z Nt Zv(N t ). (8) 

N T =0 

Due to the fact that there is a confinement condition Ylni^o^^ni — Nt for Nt in the above equation, we have 

Go (T, z) = IWo | £ exp [-PNfera] j = n ^° T3^^[Z^] - (9) 

Zq (Nt) can be obtained by noting that it is the coefficient of z Nt in the expansion of Go (T, z). Thus, we have 

Similar to the case of the developed saddle-point method for the ideal Bose gas [12,13], it is easy to get the following 
useful relations: 

n t =y: m - , (ii) 



and 
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where zo is the well-known saddle point. 

By using the above relations, the probability distribution of the system with No atoms in the condensate can be 
obtained based on the method developed in Refs. [12,13]. The normalized probability distribution function is given 
by [12,13] 



G n (N,N ) = A 

n ex P 



r N 

/ a(N,N )dN 
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where A n is a normalized constant, Nq is the most probable value of the atomic number in the condensate. If taking 
the ground-state energy of the system as the zero point of energy, a (N, N ) is determined by 
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For temperature below the critical temperature of BEC and for (3?iuj± « 1 (this condition is satisfied in the present- 
day experiments of BECs), the above equation can be approximated to be 



JVJ-JV =£ 



exp [(3e n i/ri n i) exp [-a (N, N ) /rj n i} - 1 exp [f3e n i/r] nl } - 1 
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In fact, one can get the result given by Eq. (15) through a simple physical picture. Taking the ground-state energy 
as the zero point of energy, the total energy of the system is then 



(16) 



ni 



nl 



where e e Ji can be regarded as the effective energy level of a single particle characterized by the quantum numbers n 
and I. From Eqs. (4) and (16), e e J/ takes the form 
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-nl 



£nl 



KidV + Jv^dV + Jv^dV/U 



(17) 



Thus at low temperature the system can be regarded as a fictitious non-interacting Bose system with the effective 
energy level given by e e JJ . The canonical partition function of the system reads 



^exp 



nl 



(18) 



where the prime in the summation represents the condition Y, n i^oN n i = Nt- It is easy to confirm that one can get 
the result given by Eq. (15) from this canonical partition function. 

From Eq. (15), after a straightforward calculation one obtains G n (N, N ) 
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where 
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Below the critical temperature, >> 1, the fluctuations of the condensate contributed from the collective excitations 
are given by 



(5 2 N ) = (N 2 )-(N y = ~. 



(21) 
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Generally speaking, Eqs. (20) and (21) can be used to investigate the fluctuations originated from the collective 
excitations in any dimension. Now we specify the case of 2D. For Eqs. (2) and (3), u n i and v n i are given by [15] 



U n l ~ -V n l ~ «W 2g - Xnl- 



(22) 



For a Bose gas confined in a harmonic potential V ext (r) — mojj_(x 2 + y 2 )/2, \m and e n j(= huj n i) are determined by 
the eigen equation 



(23) 



where Rj_ is the radius of the condensate. After a straightforward calculation one obtains the excitation frequency 
uj n i = ujj_ yj2n 2 + 2n\l\ + 2n + \l\, and 



(24) 



where A n i is a normalized constant determined by / |xni| 2 dV = 1. In the above equation. H n i (x) takes the form 



3=0 



where 6 



1 and the 



coefficients 



(25) 



satisfy the recurrence relation 



b j+1 /bj = (4j 2 + 4j + - 4n 2 -An- An\l\) / (4j 2 + + 8j + + 4). 

Substituting the above results into Eqs. (20) and (21), we obtain the fluctuations of particle-number in the 
condensate for the 2D Bose gas: 



NH 2 (1 - t 2 ) N^H (1 - t 2 ) N(l- t 2 ) 
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where t = T/T® with T" c ° = (N/( (2)) 1 ^ 2 huj^/k B the critical temperature corresponding an ideal Bose gas in 2D. The 
coefficients 71 , 72 and 73 are given by 
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71 ]^ (2n 2 + 2n\l\ + 2n+\l\Y 



72 



d 2 



^ (2n 2 + 2n\l\ + 2n+\l\) 



3/2 • 



73 2 (2n 2 + 2n\l\ l +2n+\l\)' (2?) 



with 



f„ (l-x 2 )x(H nl ) 2 dx , 

0nl = V L 2, ■ (28) 
J x(H n i) dx 

By a numerical calculation we obtain 71 = 0.87, 72 = 1.43, 73 = 4.37. 

From Eq. (26) we have the following two conclusions: (i) (5 2 No) is proportional to the inter-atomic interaction 
constant g. Therefore, for an ideal Bose gas there is no contribution to the condensate fluctuations due to the collective 
excitations. This is physically reasonable because the collective excitations, which are dominant at low temperature, 
originate from the interaction between atoms, (ii) The leading term for the fluctuations of the particle- number in 
the condensate are proportional to N 2 . This anomalous behaviour comes from the low-dimensional property of the 
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system. Noting that in the case of 3D Bose-condensed gas the fluctuations of the condensate due to the collective 
excitations are proportional to TV 4 / 3 [2,12]. Thus the lower the dimension of the system is, the larger the condensate 
fluctuations are. Thus at low temperature a 2D Bose gas confined in a harmonic trapping potential, which has been 
realized recently by Gorlitz et al [8] , is an ideal system for the observation of anomalous behavior of the fluctuations 
of particle-number in condensates. 

In real experiments [8,9], the Bose gases are confined in a quai-2D harmonic trap. For a quasi-2D Bose gas, the 
coupling constant is given by g w 2\/2nh 2 a s /(ml z ) [16], which is fixed by a s-wave scattering length a s and the 

1/2 

oscillator length l z — (Ti/muj z ) 1 in the z-direction, where uj z is the trap frequency in the z-direction. In this case, 
we have 



N 2 t 2 (l-t 2 ) N 3 / 2 t(l-t 2 ) N(l-t 2 ) 
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From the above equation we see that, on the one hand (5 2 Nq^ is proportional to N 2 (for large N), and on the other 

1/2 

hand it is also proportional to ui z because the factor l/l z appearing in (29). This shows that the trapping frequency 
in the z-direction plays an important role for the particle number fluctuations due to the collective excitations. Thus 
one can control the particle number fluctuations by adjusting the trapping frequency in the z-direction. Compared 
with the contribution due to the thermal atoms (see Eq. (39) in [10]), the particle number fluctuations due to the 
collective excitations are strongly dependent on the trapping frequency in the z-direction. From the relation between 
(<5 2 ./Vo) and u z , we see that the confinement of the Bose-condensed gas has the effect of increasing the fluctuations 
due to the collective excitations. This result is consistent with the role of the dimensionality in the fluctuations of the 
condensate. 

In conclusion, we have studied the particle-number fluctuations of a condensed 2D Bose gas confined in a harmonic 
trapping potential by using the probability distribution obtained through a modified saddle-point method. We have 
found that the condensate fluctuations are proportional to the inter-atomic interaction constant g and the square 
of total particle-number of the system. This anomalous behavior of the fluctuations makes it very promising to 
experimentally observe the effect of the particle- number fluctuations. The theoretical method provided here is quite 
general and can be applied to investigate the particle-number fluctuations of quasi-lD Bose-condensed gases. 
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